Nonlinear self-interaction of plane gravitational waves 
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Recently Mendonga and Cardoso [Phys. Rev. D 66, 104009 (2002)] considered nonlinear gravita- 
tional wave packets propagating in flat space-time. They concluded that the evolution equation — to 
third order in amplitude — takes a similar form to what arises in nonlinear optics. Based on this 
equation, the authors found that nonlinear gravitational waves exhibit self-phase modulation and 
high harmonic generation leading to frequency up-shifting and spectral energy dilution of the grav- 
itational wave energy. In this Brief Report we point out the fact — a possibility that seems to have 
been overlooked by Mendonga and Cardoso — that the nonlinear terms in the evolution equation 
cancels and, hence, that there is no amplitude evolution of the pulse. Finally we discuss scenarios 
where these nonlinearities may play a role. 

PACS numbers: 04.30.Nk, 42.65.-k, 95.30.Sf 
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As discussed by Mendonga and Cardoso in p| the Ein- 
stein field equations are highly nonlinear and it is a most 
interesting question whether gravitational waves share 
any similarities with other types of more well-known non- 
linear waves, for instance those occurring in nonlinear 
optics. Assuming the presence of a number of wave per- 
turbations 



h a b(n) = e ab {n)a{n) exp[iq c (n)x c 



(1) 



it was — based on the structure of the vacuum gravita- 
tional field equations — argued by Mendonga and Car- 
doso that the wave amplitude evolution equations for 
four-wave mixing processes take the form 

iq a (l)^-a(l) = «;(l)o*(2)a(3)a(4) exp[iAq b x b ], (2) 
ox a 

where (n) labels the wave perturbations, e a b(jO is the 
unit polarization tensor for wave (n), a(n) is the slowly 
varying amplitude (* stands for complex conjugate), 
q c {n) is the wave four- vector, c.c. denotes the complex 
conjugate of the preceding term, Aq b = qb{3) + q&(4) — 
5b(l) — 9fe(2) and the coupling coefficient w(l) is deter- 
mined by the Ricci tensor cubic in the wave amplitudes 
(not explicitly computed in The background space- 
time was assumed fiat. The authors then considered the 
implications of Eq. which indeed is similar to what 
arises in nonlinear optics. The analysis was restricted to 
the case when all present waves propagate in the same 
direction. In particular they focused on (i) self-phase 
modulation (when the four waves have identical polar- 
ization and wave vector) resulting in that a gravitational 
wave pulse will be continuously frequency up-shifted as 
it propagates through space and (ii) harmonic cascading 
through which gravitational wave energy is spectrally di- 
luted. 

In this Brief Report we present the results of detailed 
calculations of the amplitude evolution equation for plane 
gravitational waves propagating along x 3 = z. The per- 



turbed metric g a b is 



gab = Vab 



(,b 



ab 



■,(3) 
L ab ' 



(3) 



(2) (3) 

where n ao = diag(— 1, 1,1,1), h ab and h y ab are the second 
and third order nonlinear response to the linear wave 
perturbation . As it turns out, there are no third 
order response terms that couple back to the evolution 



equation of h^} , 



(3) 

therefore h„J is omitted from here on. 



The vacuum Einstein field equations reads 
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where the Ricci tensor has been divided into R^ b (lin 



ear in h^ b ), R { ^ b (quadratic in h}^ b and linear in h}* b ') 



? (2) 



(1) 



,(2), 



and R^ b (cubic in h^j and with terms of order h^h^). 
Having wave perturbations of the type given by Eq. 



in mind we note that in R^? , but not in R'fJ , there are 
terms that are resonant with R^j} • Thus, we identify 



? (2) 
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with the amplitude evolution equation, where it should 
be understood that only the part of R ab resonant with 



The second order 



Rab snou ld be taken into account, 
nonlinear response terms, hf b \ are related to the h^ b by 



R 



(2) 



(6) 



We are free to choose a particular gauge for h^ b and we 
choose the TT gauge such that the only nonzero compo- 

h 0- With this choice R 



nents are 



'22 



(2) 
ab 
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reads 



R 



(2) 
00 



\ [-a - d 2 t H] + \ {d t hf + hd?h, 

Ro3 = ~\dzd t H + 1 (d t h) (d z h) + hd z d t h, 



>(2) 



1 



hd z h, 



and 
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d\hf a + d z d t h^ 
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(7a) 
(7b) 
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(8a) 
(8b) 



K 2 = - [d? - dl] h 



(2) 

l a/3 — 2 L"* "zJ "a/3 

where a, /? = 1, 2 and 

i? = 4? + 4 



(ftA) - <W, (8c) 



2^4 2 3 ) + a t 2 4 2 3 ) . 



(9) 
(10) 



(2) (3) 

When computing R ab and R ab one should make use of 
the approximation d z as —dt |3|]. Applying this to R ab it 

(2) 

follows that all the components of h y ab may be set to zero 

(2) (2) 

except h\i and h 2 2 (the final result does not depend on 
whether this choice, corresponding to a freedom in gauge, 
is made or not). Making this choice Eq. © becomes 



d 2 H 



(d t h) 2 + 2d 2 (h) 2 



(3) 

The nonzero components of R ab are 

Mi ) = \{d z h){d z H)- ~(d t h)(d t H), 



(11) 

(12a) 
(12b) 

(3) 

and with the approximation d z « — d t we see that R ab = 
0. The amplitude evolution equation JSJ thus reduces to 



R22 = --(d z h){d z H) + -{d t h){d t H), 



i 



i 



(d*-d?)h = 0, 



(13) 



and we conclude that — up to third order in amplitude — 
plane gravitational waves propagate with no change in 
the amplitude modulation, i.e., w(l) in Eq. @ is zero. 
To some extent this may seem counterintuitive, since a 



gravitational wave-packet propagating in flat space-time 
transports energy and momentum as described by the 
pseudo energy- momentum tensor which is of order h 2 . 
This produces a local curvature of space-time of the same 
order, and one could expect a coupling to the original 
wave perturbation, leading to cubic nonlinearities in the 
amplitude evolution equation. However, this picture is 
only partially true. The pseudo energy-momentum ten- 
sor does produce curvature quadratic in the amplitude h 
(essentially described by Eq. fTT)) '). and this couples back 
to the original perturbation leading to cubic nonlinear- 
ities in the Riemann tensor. But as follows from Eq. 
(I12all and Eq. (|12b|) these terms cancels when computing 
the Ricci tensor and thus does not enter the amplitude 
evolution equation. 

It is instructive to compare this result to existing ex- 
act plane wave solutions of the vacuum Einstein field. 
One particular exact plane wave solution is given by the 
metric (see, e.g., Ref. £|) 

ds 2 = -dt 2 + a(£) 2 dx 2 + b{£) 2 dy 2 + dz 2 , (14) 

where a and b satisfy the equation 

6<9f a + adp = 0, (15) 

where £ = z — t. Focusing on weak waves, we take 
a = (1 + h+ ^H)i and b = (1 - h+ ^H)i where 1 > h 
and H ~ h 2 . Expanding Eq. (|15|) up to third order in 
amplitude one find precisely the relation (| 1 1 11 between H 
and h and that the third order terms cancels identically. 
The fact that the metric components of the gravitational 
wave perturbations in Eq. (with the gauge choices 
made above) coincide with those of the exact wave so- 
lution when expanded in powers of amplitude suggests 
that Eq. I|13|) will remain even if one relaxes the approx- 
imations made in this Brief Report. 

It should be pointed out, however, that gravitational 
four-wave mixing processes are not ruled out in gen- 
eral. Firstly, we have here only considered parallel prop- 
agation, and it is well-known that anti-parallel gravita- 
tional waves can interact (see e.g. Ref. y| and references 
therein). Secondly, the presence of background curva- 
ture and matter modifies the properties of gravitational 
waves, e.g., giving rise to wave dispersion and new types 
of nonlinearities. To what extent this increases the pos- 
sibilities of nonlinear self interaction is an open question. 
The equations presented here provides a starting point 
for future investigations of these possibilities. 
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